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MSTRACT 

During the initial space-time evolution of an electron beam injected 
into the corona, the strong beam-plasma interaction occurs at the head 
of the beam, leading to the amplification of a quasi-monochromatic large- 
amplitude plasma wave that stabilizes by trapping the beam particles. 
Oscillation of the trapped particles in the wave troughs amplifies 
sideband electrostatic waves. The sidebands and the main wave subsequently 
decay to observable transverse electromagnetic waves through the parametric 
decay instability. This process gives rise to the elementary striation 
bursts. Owing to velocity dispersion in the beam and the density gradient 
of the corona, the entire process may repeat at a finite number of discrete 
plasma levels, producing chains of elementary bursts. All the properties 
of the type Illb bursts are accounted for in the context of the theory. 
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I. INTRODUCTION 

The so-called type Illb solar radio burst was defined by de la Noe 
and Boischot (1972) as a chain of several elementary bursts which appear 
in dynamic spectra as either single, double, or triple narrow-banded 
striations (Figure 1) . The spectral features of these striations were 
described extensively by Ellis and McCulloch (1966, 1967) and de la Noe 
(1975) . V7e shall describe here the salient characteristics of the 
bursts, which will be theoretically explained in this paper. 

A. Properties of Elementary Bursts 

1. Individual striations exhibit a very narrov? bandwidth Af: 
the average is approximately 100 kHz, but Af may vary from the frequency 
resolution of the receivers ('v- 15 kHz) to a few hundreds of kHz. 

2. In double- or triple-striation groups, the frequency separation 
(Sf , between the edges of the individual striations, is also approximately 
100 kHz on the average, with less-dispersed limits than for Af. 

3. The duration at a given frequency varied from the time 
resolution of the receivers (2^ 20 ms) to a few s; the average duration is 
about 1 s, a factor of 4.5 lower than the duration of normal type III 
bursts (Aubier and Boischot 1972, de la Noe 1975). 

4. These short, quasi-monochromatic bursts are observed at 
frequencies between 100 MHz and 10 MHz. No such bursts have ever been 
reported at higher frequencies. A search by one of the authors (J.L.N^) 
to detect them at lower frequencies among satellite observations by IMP-6 
and RAE-2 was negative , possibly because of the frequency and time 
resolutions of these experiments. The occurrence of elementary bursts 
was found to peak around 30 MHz (de la Noe and Boischot 1972) , at least 
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for the studied period July-August 1970. 

5 . Circular polarization is observed for the striations , but 
it occurs in a totally random way from one to another elementary burst in 
a chain. However, when doublet or triplet elementary bursts are found 
polarized, all the striations have the same sense of circular polarization. 

6. In a few cases 20 %) , elementary bursts were found to 
have a non-zero drift rate towards lovjer frequencies.. When it occurs, 
this frequency drift affects most striations occurring in a short period 
of time. The drift rate may reach values up to -150 kHz/s, but in most 
cases is null or not measurable. 

B. Properties of Type III b Bursts 

Besides the intrinsic properties of the elementary bursts, 
described above, the properties of chains of elementary bursts lead to 
interesting clues to understanding the phenomenon. 

1- The leading edge of the spectral envelope of type Illb 
bursts is similar to that of type III bursts. This suggests not only 
their appellation, but also that the type Illb bursts are produced by 
the same kind of exciter as type III bursts: viz., electron streams with 
average velocity 2; c/3 (Lin, Evans, and Fainberg 1973). 

2. de la Noe and Boischot (1972) found that 30% of t 3 /pe Illb 
bursts were precursors to type III bursts, strengthening the hypothesis 
of a common agent for both types of bursts, de la Noe (1974) showed 
that this relationship, which could be owing to pure association, was 
real mainly when the sources of the bursts are located near the limb. 
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although location of sources on the central part of the disk does not 
preclude such association. Type Illb bursts seem to be emitted at a large 
angle to the local magnetic field, rather than along the field. 

3. A recent study of type III bursts by Rosenberg (1974) led 
that author to conclude that 95% of all type III bursts occur at 2 
where oig is the local plasma frequency. When a type Illfo burst is seen 
as a precursor to a type III, the measured frequency ratio 

at a given time is somewhat dispersed. It suggests, however, that emission 
of the type Illb burst takes place close to Wg, while the type III is at 
2o)^. This assvimption is strengthened by the recent observation of both 
fundamental and harmonic emissions of a type U burst, in which the fundamental 
is striated while the harmonic emission is smooth (Stewert 1975) . 

4. Positional measurements of precursor type Illb and following 
type III bursts showed that both emissions are located at the same place 

in the corona (de la Noe and Gergely 1975) . This strengthens the hypothesis 
that both bursts are excited by the same electron stream with two different 
mechanisms: the first produces the type Illb burst, takes place at the 
front of the beam and requires a shorter time to generate the emission; 
the second one which generates the type III burst is a process like that 
studied by Papadopoulos , Goldstein, and Smith (1974). 

C. Outline of the Theory 

We assume that the particles are injected from a localized 
acceleration region in the lower corona, forming a beam which is guided 
by the magnetic field. The velocity dispersion in the beam causes a time- 
dependent evolution of the distribution function at any point in space. 

For some time after the beginning of the injection, no particles will have 
arrived at, say, the point x. Then the fastest particles in the beam 
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arrive at x, followed by successively slower particles. For some time during 
this evolution, the quantity Vj^/U is very small, where is the thermal 
spread that characterizes the velocity range over which the beam distribution 
rises to its instantaneous peak at the velocity U. In the theory that follows, 
we shall assume 

u 

where r) H nj^/n^ is the instantaneous ratio of the beam density to that 
of the ambient plasma. Although the condition (1-1) may seem stringent, it 
is quite easy to construct models ox the beam evolution in which (1.1) is 
satisfied. To preserve the continuity of the present discussion, we defer 
to Appendix A a fuller examination of the evolution of the beam. Here 
and in what follows we shall simply assume that inequality (1.1) is satisfied 
at some time in the early evolution of the' beam at a given point. In 
Appendix A we show that this assumption appears justified under a wide range 
of plausible conditions for the injection of the beam into the corona. 

We may now discuss the salient features of the theory that follows. 
When a beam satisfying inequality (1.1) propagates through a background 
plasma, quasi-monochromatic plasma waves, with phase velocities just below the 
beam velocity U, are unstable with exponential growth rate y 

We refer to this regime as the strong beam-plasma interaction, in contrast to 
the so-called "bump-on-tail" instability in vhich inequality Cl.l) is reversed 
and a broad band of waves grows at approximately the rate hiOg,. In the 
strong-interaction regime, the growth of the beam-mode waves is saturated 
by trapping of the beam particles . The oscillation of the trappecl particles 
in the wave troughs gives rise to unstable si<leband wayea, geparated in 
frequency from the main wave by approximately the bounce frequency of the 
particles. The main wave and the sidebands^ which are all electrostatic 
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oscillations, subsequently decay into transverse waves, at approximately 
the plasma frequency, and ion-acoustic waves. In addition, as the side- 
bands grow to amplitudes comparable to that of the main wave, the particles 
experience a demodulation of the wave phases; i.e., the autocorrelation 
time becomes of the order of the bounce period. The particles are then 
no longer trapped. After a certain distance, however, owing to the 
velocity dispersion of the beam and the density gradient of the corona, 
the beam will again satisfy condition (1«1) and the entire process may repeat. 

V7e now discuss each facet of the theory in turn. 
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II. THE BEAtl-PLASMA INTERJ.CTION 

The interaction of a cold low-density beam with a plasma has been 
investigated by O'Neil and Malmberg (1968) , Drummond ad . (1970), and 
O'Neil, Winfrey, and Malmberg (1971) . We shall consider an electron beam 
in,jc .ted into a neutral background plasma. This system is not charge- 
neutral, and a return current flows with velocity gU. As was shown by 
Melrose (1974), the effects of the return current are negligible for 
conditions appropriate to type III bursts , and so we consider the ions to 
be a fixed neutralizing background. As usual, we consider a one-dimensional 
interaction. Then the dispersion relation for electrostatic waves is 


1 


f dv k 8F„/9v 

-F- j 



k 9Fb/3v 
kv-w 


0 , 


( 2 . 1 ) 


where o)g is the electron plasma frequency, and Fg, Fj^ are the distribution 
functions of the ambient and beam electrons, respectively. 


Following O'Neil and Malmberg (1968) , we take the beam distribution 
to be Lorentzian: 


Fb (V) 



1 

(v-U)^ + 



( 2 . 2 ) 


where measures the velocity spread. The assumption of the form (2.2) for 
Fj^ is a matter of convenience; O'Neil and Malmberg showed that results similar 
to those below are obtained by taking Fj^ to be Maxwellian. The principal 
advantage of the form (2.2) is that it enables one to evaluate the resonance 
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integral for the beam in Eq. (2. l)easily*; 


^ 9Fb/9v = (w - ku + ik Vj^) 
I kv-u3 


-2 


(2.3) 


Vie denote the plasma dielectric 


J 


e^(k,£i3) = 1 - j" dv k dF^/dv . 


kv-u) 


Then, assuming 


exp {-U^/2V 2) << ri. 


(2.4) 


3 . 

where Vq = (T^/m) is the electron thermal velocity, e(k,w) has an 

o 

approximate zero at the point 


= u)e (1 + 3 k^Vg^)^, 


^0 “ Wq/U. 


(2.5) 


Inequality (2.4) Is easily satisfied under conditions of interest to us, in 
which we assume 5 ^ U/Vg ^ 20 and 10“® ^ n 10“^. Then, using (2*3), 
the dispersion relation (2.1) may be expanded about (w^,k^) to give (O'Neil 
et al. 1971) 



where m' = u - oJq and k' = k - kQ. In the present context, the waves 


*We note that Eq. (2.2) implies an infinite energy density for the beam. This 
will be of no concern to us, since we shall not take moments of the distribution 
(2.2) but instead merely take advantage of its analytic properties. Note that 
(2.2) is a representation of the Dirac delta-function in the limit 0. 
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vary as exp i Ckx - wt) , so that Im u* > 0 implies instability. The 
integral in Eq Ck,(D) is to be performed on the usual Landau contour. 
Then we have 


Re e Ck,w) = 1 - 


v2 1 2 _ 


3 X 


De 


where Xjjg is the electron Debye length, and so 


0)0 /'SEoi 


\ 9(0 ■ 


k 3e 
o I 0 

9k 


- 2(1 + 3C2) , 




(2.7) 


(t»o,ko) 


where C = Vg/U. Differentiating (2.6) with respect to k’ , using (2.7), and 
maximizing = Im w' , one finds that maximum growth occurs at k' = 0, 
where the real frequency w and growth rate are given by 


0) = Wq (1 - , 

Y = (3^/^/2^/3) ;!^l/3 


( 2 . 8 ) 


where n = n (Wg^/Uo^) [(%/2) (9eo/9^^)w ,kj] Note that , given 


by Eq. (2.5), contains the Bohm-Gross dispersion for a warm plasma. In this 
1/3 

limit (V./U < ri '^ ) , the instability is reactive. The group velocity V 


Vg - 

dk 


a u + i Vg2 
^ ^ u 


(2.9) 


is given by 
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Furthermore, the half-width in k-space of is given approximately by 


Ak - 1 


8 ^ 


3k" 


Ib. 


Y=Y max 


= 3 

J5/6 


-vl/3 , 
n k 


( 2 . 10 ) 


We consider that the unstable waves grow out of the small thermal 
noise level present before the beam enters the plasma. Then many growth 
times will pass before nonlinear effects become important. After N 
e-foldings (t = N Yj^"^) r the half-width of the unstable wave spectrum will 
be 

6k = (In 2/n)^ Ak; 

i.e. , 

|E(kjj ± 5k) 1^ = 1 |E(ko) 1^. 

2 

The waves stabilize by trapping particles when their amplitudes become 
sufficiently large to alter significantly the reactive orbits of the 
particles (O'Neil 1965, O'Neil et al. 1971). This process occurs only 
over the last few e-foldings. The question then becomes whether the 
spectral width 6k is sufficient to destroy the autocorrelation of the 
spectruiTi in a time y.”^* phase difference seen by an electron of 

D 

velocity v in this time interval is given by 

A4, = I V - Vg I Yb”^ 5k, 

and for autocorrelation to be maintained we require A(p«l. For the plasma 
electrons this condition becomes 


N** >> [scin 


( 2 . 11 ) 
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and for the beam electrons 


N 


h >> [2Cln 


( 2 . 12 ) 


Subsequent to trapping, the beam particles exchange energy with the 
wave; the exchange of energy occurs with the bounce frequency 


“b “ (ekoEg/m)^, 


(2.13) 


where is the saturation amplitude of the wave. Violation of inequality 
(2.11) causes nonresonant velocity diffusion of the background electrons, 
while violation of (2.12) causes resonant velocity diffusion of the beam 
electrons. For our purposes only the latter is of direct concern. We 
shall see below that under conditions of interest to us, the sides of (2.12) 
are in the ratio of order 5:1. Thus, it might appear that N~^ is only 
marginally valid as a small parameter. Nevertheless, several recent experiments 
and theoretical studies indicate that the single-wave theory is valid 
despite the implications of inequality (2.12). in order not to interrupt 
our exposition here, we present a discussion of this point in Appendix B. 

In what follows, we assume that the single-wave theory is valid 

for times of the order of several bounce periods after saturation of the main 

wave , 

At saturation, the time-averaged energy density of the trapping wave 
(averaged over one period) is given by 


Eg^ ^ 2"^/^ Cn. m UV2) - 2""V3 ^4/3 jj2 


16 -rr 


V, 




(2.14) 


where Tg is the temperature Cih eVl of the background plasma. 
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III. THE SIDEBAND INSTABILITY 

The motion of the particles trapped in the large-amplitude-wave 
troughs is significantly distorted from the trajectories given by linear 
analysis. In the instantaneous frame moving with the phase velocity of 
the trapping wave, the trapped particles execute closed orbits in phase 
space (Pig. 2) , whereas untrapped particles move unidirectionally over 
the wave with a sinusoidal modulation of their velocities relative to 
the wave. The amplitude of this modulation is 

Vos ® ^s/i" (3.1) 

The distortion of the orbits of the trapped particles leads to 

the coupling of the main wave to sideband waves at frequencies coq ± Wg, 

where the bounce frequency a)_ is given by Eq. (2.13) . The parametric sideband 

s 

instability has been examined by Kruer, Dawson, and Sudan (1969) , 

Goldman and Berk (1971) , and Wong (1972) . Other recent treatments have been 
given by Bud'ko, Karpman, and Shklyar (1971) , and by Brinca (1972) . In this 
section we shall present only the salient points to be used in the remainder of 
this work; additional remarks about the sideband instability are given in 
Appendix B. 

The initial observation of such sidebands was reported by Wharton, 
Malmberg, and O'Neil (1968). In the experiment by Wharton et al., a large- 
amplitude wave was launched into a plasma. In such a situation, the trapped 
electrons separate generally into two classes: "resonant" particles 
which are trapped near the bottom of the wave trough, and "nonresonant" 
particles which oscillate on trajectories near the separatrix of Figure 2. 
Generally speaking, the resonant particles execute closed trajectories in 



phase space with the bounce frequency ciJ_, while the nonresonant particles 
take an arbitrarily long time to execute a closed orbit. The treatments 
of Kruer et al. (1969) and Wong (1972) are appropriate to this situation. 

In our case, however, we are dealing not with a large-'amplitude 
wave launched into the plasma, but rather with a wave which grows from 
thermal noise on a beam. In this situation, nearly all the particles 
are resonant. The theoretical treatment appropriate to this case is 
that of Goldman and Berk (1971). 

Let us denote the frequency and wavenumber of the trapping wave 
by and , respectively. Then the principal sideband waves are at 
frequencies co = Wq + Aw and wavenumbers k = k^ + Ak which are given by 
the dispersion relation of Goldman and Berk (1971)*: 



where e = n under the assumption that all of the beam particles 

are trapped. From Eqs. ( 2 . 13 ) and (2-14) , we find 

e = Ti (n Eg^/4ir m = 2. 


*Here, Ak should not be confused with the Ak of Eq. (2.10). 
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In Figure 3 (a) we have plotted the maximum sideband growth rate , 
and the wavenxunber shift and frequency shift at which the maximum growth 
rate occurs, as a function of ^ = Ve/U; all quantities are normalized to 
(i)g. We see that all of these quantities are only weakly dependent on 
i in the range 0.1 ^ ^ 0.2, while they approach asymptotic values for 

5 ^ 0.1. (For a coronal temperature of 2 x 10 K, a value of ? = 0.1 
corresponds to a beam energy m 0^2 = 20 kev, while ^ = 0.2 implies 
an energy =5 kev.) Furthermore, the results are found to be only weakly 
dependent on e, in agreement with the results of Goldman and Berk. In 
addition, denoting AtJ = + i Ygb' solutions obey the symmetries 


Ao)j. (-Ak) = - Ao);^ (Ak) , 


■fsb * ^sb • 


(3»3) 


The growth rate Ygj^ drops off sharply for |Ak| > ^ 

tapers off more slowly for |Ak| < |Ak(Y ) |; the width in frequency Aw 
between half-maxima is approximately Wg. Analogously to the discussion 
of^II regarding the narrowing of the spectrum as the waves grow, the 
frequency width of the sidebands becomes ^.2 w as the most unstable 
wave becomes comparable in amplitude to the main wave. 

In Figure 3(b) we have replotted the results of Figure 3 Cal normalize<i 
to Wg for various values of n? from Eq. (2 -141, the normalization is given by 
“ 2-V12 ,1/3 
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IV. THE PARAMETRIC DECAY INSTABILITY 

Large-amplitude, quasi -monochromatic plasma w^ves are unstable to 
fixed-phase wave-coupling interactions known as parametric instabilities. 

In general, such instabilities may occur when there exists an approximate 
frequency-matching relation between the large-amplitude "pump" wave and 
two or more normal modes of the plasma. In contrast to weak-turbulence , 
random-phase interactions, however, the unstable modes may have dispersion 
properties different from those of the normal-mode waves; i.e., they may 
be eigenmodes only of the nonlinear plasma dielectric in the presence 
of the pump wave. A case in point is the oscillating— two-stream instability 
(Nishikawa 1968a/b; Sanmartin 1970) , in which the pump is a dipole Ch=0) 
plasma wave at the plasma frequency, the normal modes are a plasma wave 
with frequency higher than the pump wave Ch>0) and an ion acoustic wave, 
and the unstable modes are frequency-shifted ion-acoustic and plasma waves 
with|k|>0. This instability was applied to the problem of the stabiliza,tion 
of the type III exciter by Rapadopoulos et al. C19V4) . 

Parametric instabilities are also characterized by the existence of 
a threshold amplitude which must be exceeded by the pump wave in order 
that instability may occur. In the present context, the only instability 
for which the threshold is exceeded is the decay of the Langmuir pump wave 
into a transverse electromagnetic (TE) wave and an ion-acoustic oscillation. 
This process was investigated by Lashmore-Davies (1974a, b). 

Consider a Langmuir pump wave, a TE wave, and an ion-acoustic wave 
with frequencies and Wg, respectively, and wave vectors kj^, ^d 

it , respectively, such that . 
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L T ? ' 


0) - 03 +0) 

L T S 


(4.1) 


Denoting 6 S co^ - oj^, we look for unst^le modes with, complex eigenfreguencies 
0 ). Further denoting = (o + fi , the dispersion relation for a traveling 
Langmuir-wave pump was shown by- Lashmoret^Davies C1974a/b) to be given by 


+ i (Si" 


03c 


+ 2 i Y “• -4 K «e Of 


(4.2) 


where Y„fY„ are the damping rates for the TE and ion-acoustic C.SJ waves, 
respectively, and where 


K 1 


U3 


(4.3) 


4 8tt n^T^ 


The dispersion relation (4.2) also describes the up-conversion of a 
plasma wave to a TE wave with higher frequency. In this case the matching 
relations are ^T ~ ^ 

thresholds for the two processes are (LashmorerDavies 1974al 

, 5<0 ; 


= * Yg |s|/“s 


K = 2 6 03„ 

c s 


, 6>0, 


(4.4) 


In Eq. (4.4), the condition [ 5 [ > C.2 Y^Yg) must also be satisfied* For 
Y_ - 03_, the instability threshold for both 6 ^ Q is K = 4 YeYrnf independent 
of 6. Therefore, we must estimate the relative magnitude of Y,p/“g and the 
value of K in both the main wave and sidebands. 

To estimate K we assume that the sidebands grow to approximately 
the amplitude of the main wave. This assumption is suggested by the 
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observation that when multiple- striation elementary bursts appear, all 
elements in the burst group appear to be of roughly equal intensity. 
Furthermore, owing to the presence of the sideband waves, the particles 
experience wave phases of several times 2ir during one bounce period; i.e. , 

6co = l^sb “ ^0 I ^ Therefore the particle motion becomes stochastic, 

rather than coherently modulated as in the trapping wave. Without examining 
the exact particle dynamics, it is reasonable to assume that this demodulation 
becomes important when all the waves are of approximately equal amplitude. 
Therefore, we estimate K for the sidebands as well as the main waves by 
substituting Eq. (2.14) into Eq. (4.3). The resulting value is 



where the subscript "sat" indicates that the saturation amplitude of the 
main wave has been used. 

The value of to be used in determining whether or not exceeds 

is the enhanced damping decrement determined by the presence of the 
ion-acoustic turbulence that is produced by the decay instability itself. 

If, however, the TE wave has a frequency which lies within the frequency 
range of the driver wave packet — either in the main wave or the side- 
bands — then y^ should be replaced by 6o3, which is an effective collision 
frequency. In Appendix C we estimate y^ and show that with such considerations 
the threshold values (4.4) are exceeded by K only for 6<0. In Figure 4.4 

S3.U 

we plot K and K ^./K vs. 5 for 10“® < n < 10~^. We see that the 

decay threshold is typically exceeded by factors of order 10-1000. It may 
therefore appear at first glance that K is not the correct pump intensity 

S 3.4- 
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to use, because the decay instability will occur at smaller levels as the 

p\amp is growing. We shall see below, however, that even at the saturated 

amplitudes the decay rate is much smaller than y^^ and Ygjj* Thus the 

first three phases of the elementary burst process — growth and saturation 

of the main wave, growth of the sidebands, and decay of all the waves — 

occur on well-separated time scales. 

We denote the wavenumber of the TE wave by k^, and define = ck^u^ 

Because a « 1, the index of refraction 
T 


ck a a^(l - a^V2) ~ a^. 
0) 


Now consider the decay of a pump wave of phase velocity V, , and let 




V s V /V , . Then the frequency mismatch 6 is given by 
e (p 


£ = (1 + _ (1 + 3 v 2 )’*. 


(4.6) 


and 


1 , 3 _ 

^ = y^(V “ Ci^ V^/C) - V. 


(4.7) 


Prom (4.6) , the range of 6 is 


0 < < (1 + 3v2)^ - 1. 




As expected from the frequency matching condition + oig, the growth 

rate of the decay instability is a maximum for |5| = Ug. 
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In Figure 5 we show (dashed curves) the maximum growth rate, calculated 

from the dispersion relation (4.2), as a function of V for various vtlues 

-5 -3 

of n. This rate is typically in the range 10 < y,/o) < 10 , and the 

functional form of the dashed curves in Figure 5 is 


^d -2 ^sat 

— = 1.1 X 10 5 

(D U) 


( 4 . 8 ) 


e s 

The dispersion rela.tion (4.2), however, assiames a monochromatic pximp 

spectrum, an approximation valid only as long as the growth rate is 

greater than the frequency width 6o) of the pump. We have seen above that 

6u ~ 0.2 u , and in general the condition of a monochromatic pump for the 

decay insteibility is invalid. In this case, we must modify the growth rate by the 

prescription (Valeo and Oberman 1973) 


Y^^/5co, 


and so we find 


M .. ,..-4 ^ sat 

- — 5 X 10 - -— 7 - — jr 

(b 1/3 u)^ 

e n s 


(4.9) 


The corrected value of y^ given by (4.9) is shown in the solid curves of 
Figure 5. From Eqs. (4,6) and (4.7) we see that the index of refraction 
of the TE wave for which y^ is maximum (Figure 6) is independent of 
n and tends sharply to zero for V < 0.05. Therefore, we shall see below 
that only the range V > 0.05 is of interest to us. The width 
of the decay spectruip, however, varies with n, ar»d therefore a finite 
bandwidth in is produced. In accord with our earlier convention, we 
define the width of the decay spectrum as that part of the spectrum for 
which y^ > .5 y^ max* bandwidth of the resulting TE spectrum is 

A(o _ Aa^ 

T 
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In Figure 7 we -plot this quantity. In the next section, we discuss these 
results in the context of the type Illb observations. 
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V. APPLICATION TO TYPE Illb BURSTS 

As Noted in<^I, the characteristics of elementary bursts vary 
widely about certain averages. VJe may understand this variability by 
noting that the growth rates and spectra of the various instabilities 
involved in the theory are somewhat sensitive to the instantaneous density 
and velocity of the beam (or, equivalently, the phase velocity of the 
decaying sidebands) . These features are in turn determined by the overall 
intensity and shape of the beam v^iiocity profile that is injected into 
the corona (see Appendix A), and these parameters may in principle vary 
widely among events . The corona itself is also variable in time , with 
regard to such features as its temperature and density profiles. Therefore, 
we can realistically ask of a theory for the type Illb bursts only that 
it offer a framework within which the average phenomenological features 
of the bursts, and their range of variability, may be accounted for. In 
this section we show that the present theory meets this test. We shall 
focus first on the properties of the elementary bursts, and then consider 
how chains of such bursts may occur. 

A. Properties of Elementary Bursts 

1 . Frequency Separation of the Striations 

When two or three striations occur in an elementary burst, 
the frequency separation between them should be approximately the bounce 
frequency w^, according to the theory of the sideband instability (<^III) . 

At 30 MHz (the frequency of peak occurrence for type Illb bursts), the 

—3 

average frequency separation (Sf is around 100 kHz. Thus §f/f 3 x 10 f 
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equating 

f tOg 

^6 

gives an average value of 3 x 10 for n* Allowing for a factor of about 

3 on either side of the average 6f (which is not widely dispetsed) leads 

. -9 -6 

us to consider the density range 10 < n < 10 ; the limits are extremes 

—8 —7 

and the observations indicate the narrow range 10 ri 10 as most 

—8 -"7 

typical. Thus we have presented numerical results for n = 10 and 10 

* 

That this is a reasonable range of n can be inferred from considerations 

of both the f\indamental condition of the theory, given by inequality (1-1) , 

and the contrasting requirements necessairy for the development of the 

normal type III burst. Turning first to the latter of these considerations, 

we note that if (1-1) is violated, the growth rate of the weak beam-plasma 

instability is - n Wg. The collisional damping rate is where 

3 8 

ue typically of order (1-3) x 10 in the lower corona. Therefore, 

-»8 *7 

a beam of instantaneous density 10 h 10“ ' , and in which condition 
(1-1) is violated, will be too weak to drive the linear beam-plasma 
instability. Conversely, (1-1) is most likely to be satisfied precisely 
at those times in the evolution when only the fastest particles have 
reached any particular point, and these particles constitute only a small 
fraction of the injected distribution. 

2 . Bandwidth and Frequency Range of Striations 

As can be seen from Eq. (4.7), the refractive index of the 
TE wave is a function of the phase velocity of the decaying pump wave. 

In 41 we noted that the typical bandwidth Af of a striation was of order 100 
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kHz, ranging from <15 kHz to a few hundreds of kHz. At 30 MHz, this is 
a range of 

5 X lO’^ Af/f < 10"^, 

with an average Af/f of about 3 x lo”^. Prom Figure 7, we see that this 
implies 


0.06 < V < 0.20 , 

'Vi 'V, 


(5.1) 


where V s V^/V^, For a coronal temperature of 2 x 10 K, this corresponds to an 
energy range of 5 < < 80 kev, if the phase velocity is assumed to be equal 

to the beeim velocity U. It is doubtful that electrons with energies as 
low as 5 kev actually play a role in type Illb bursts, because under most 
conditions (1.1) is probably violated by the time the beam evolves down 
to such low energies. Figure 7 indicates, however, that whether or not V 


<!> 


is equated to U, we need to consider V, in the range (5.1) 

9 


The role of V is fundamentally to determine the frequency 
of the pump via the Bohm-Gross dispersion relation. Therefore the sidebands, 
which are shifted from the trapping wave by about twice the bounce frequency, 
are described for purposes of the decay instability by effective phase 
velocities different from that of the central trapping wave. Denoting by 
V^, and V the ratios of to the phase velocities of the main wave, upper 
and lower sidebands, respectively, we have 


= 1 


(1 3 V + Aur 

/3(L ° "ST 



(5.2) 


where AWj. is given as a function of Vq and H by Figure 3a. In Figure 8 
we show corresponding to Vq for various n . Note parficularly- that at 
low values of Vq, V_ may be < 0.05 ^ As noted above, at V the refractive 

index of the daughter TEu waves decreases rapidly (Figure 6) . ■ This is of 
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great significance when we consider the question of why elementary bursts 
most commonly consist of one striation, less frequently of a pair of 
striations, and only very occasionally of a triplet group. 

3 . Single- vs. Multiple-Striation Bursts 

To understand this, let us consider the attenuation of the 
transverse waves by collisional absorption in the corona. For frequencies 
near the plasma frequency, the absorption per unit length is given by 



Vq 


(5.3) 


Because the wave is propagating through a medium of decreasing density, 
the optical depth t may be approximated by 

T K L, (5.4) 

where L = ngjdng/dr( is the scale length of the coronal density gradient. 
Thus 

T - '^o/“T' 

where Tq = Vq L/c. In Figure 9 we show values of L and derived from 
three coronal density models which represent a wide range of coronal 
conditions, ranging from the quiet corona at solar maximum (Leblanc, Leroy, 
and Pecantet 1973) to Newkirk's active streamer model (Newkirk 1967)^ 
to an intermediate state (Riddle 1974) . We see that Tq ranges from values 
of 0.1 at 30 MHz to ^ 0.4 at 100 MHz. On the other hand, we see from 
Figure 6 that pump waves in the range 0.05 ^ V ^ 0.15 decay to TEm waves 
with 0i07 < 0.20. Therefore, at any frequency, the higher the frequency 

of the decay p\amp, the higher the index of refraction of the daughter 
TE wave and the lower the attenuation factor exp(-TQ/arj,) . As a specific 
illustration, we combine in Figure 10 information from Figures 6, 8, 9 to 
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show the optical depth of the three possible elements in a triplet, at 
30 MHz and 100 MHz, using Riddle's density model and an assumed density ratio 

-7 

n = 10 . The main wave phase velocity can be equated to U, so that 

Vq = 5. We note that at 100 MHz the waves suffer severe attenuation at 

all V„. At 30 MHz T becomes <1 for the lower sideband only at V =0.09, 

o o 

corresponding to beam energies of =25 kev. 

+ o ~ 

At all frequencies and at all beam velocities , t < x < t 
( where the superscripts o, +, - refer to the main wave and the upper and 

lower sidebands, respectively) . Depending on the frequency, however, one 

or more of the waves may suffer severe attenuation by exp(-x) between 

the source and the observer. For example, note that although both the 

main-wave and the upper-sideband daughters suffer little attenuation at 

30 MHz, a similar calculation at 60 MHz would lead to x°<l only for 

Vq= ? > 0.10, while the corresponding value of x'*' is 0.8. 

The optical depth, especially defined only approximately 

as in Eq. (5.4) , is admittedly a rough measure, and these arguments are meant 

to be suggestive rather than definitive. Nevertheless, because under 

most conditions we would expect condition (1-1) to be satisfied principally 

in a narrow range of probably bounded above by ? < 0.10, the consideration 

of the attenuation of the various members in the elementary bursts offers 

a very plausible explanation of why single striations — corresponding to 

the upper sideband only — are seen about 80% of the time, while doublets — 

decay products of the upper sideband plus main wave — are seen about 20% of 

the time, and triplets, which require unusual conditions to be met by the 

beam, are rarely seen. It also explains why type Illb bursts are not observed 

at frequencies greater than about 100 MHz; at such frequencies, is 
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relatively large, and severe attenuation of the emission would be expected 
to occur. CAn explanation is given below for why bursts are not observed 
at lower frequencies.) 

We have stressed the variability in the bursts owing to 
the variability of the beam evolution, because of the fundamental importance 
of condition Cl-1) and because the dependence of the instability parameters 
on 5 = Vg/U is easily seen. We note, however, that variations in the 
corona, such as temperature variations, for example, may also be important 
in determining the local parameter regime. An indication that this may be 
the case is the fact that despite the arguments above, which suggest that 
multiple-striation bursts ought to occur more commonly with decreasing 
frequency, we have not found a strong correlation between frequency and the 
niamber of elements in a burst group. This could be due, for example, to 
the decrease in temperature with decreasing plasma frequency, which would 
tend to keep C = Vg/U at low values and suppress the lower sideband. 

4- Duration of the Elementary Bursts 

Let us turn to the topic of the observed duration of the 
burst at a given frequency. We assume for the moment that this duration 
is equal to the duration of the decay- instability phase, which we denote 
by t. In I we noted that the observed duration varies from less than 
20 ms to a few seconds, but that its average value is about Is. As noted 
in ^ IV, we may write the maximum growth rate of the decay instability as 

» 5 X 10 K “e/ (5.5) 




where K = K/u Then the equation for the normalized TE energy density T 
s 


IS 


,'V« 

S - V, (K)T, 


(5.6) 
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while 

^ % 

^ / {5.7) 

dt dt 

where in (5.7) we have neglected the energy of the ion-acoustic waves 
compared to that of TE waves. Analogously to K, T xs defxned by 

T S ^ 

4t0s^ ( 

'll 'X; , , 

so that K/T gives directly the ratio of energy densities in the pump 
and daughter waves. In Eq. (5.6) have also neglected the convective 

transport of TE waves out of the source region, because its rate is small 
compared to Equation (5.7) has the solution 

'll '1/ 'll 'll 

K(t) = Kq + Tq - T(t), 



where subscript zero denotes initial value at the beginning of the decay . 




Substituting this into Eq. (5.6) ,we may solve for T(t) in quadrature; 


3 

U t = 2 X 10 — rj— 

e 


f\j '\, 

K + T 
0 0 


tn 


r 

T/T 

0 . 


a. 

1 + T 

0 

- T 

Sf 

0 

'V. 

K 

0 


(5.8) 


Then denoting by t the time at which T = K 


o 
1/3 


, we find (using ^ « K ) 


00 t 4 X 10' 
e 


TZhT 
^ K 


In 


r ^ . 
0 


'\i 

T 


(5.9) 


0-' 


Assuming that, the initial noise level in TE modes is of the same order 
of magnitude as that in electrostatic modes, we may employ the arguments 
of Appendix C to estimate 

^ . 10® , 

5'o 5’ 
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where we replace by unity because the dependence on this factor is 

only logarithmic. Furthermore, taking the initial pump level to be given 

by K we have K * where u == 1/1836 is the electron-to-proton 

sat 0 

mass ratio. Defining e by C = 0.1 e, with e of order unity, we finally 

obtain w t = c^/5r\, so that 

e 

- ^ 3 X 10~^ (s) ^ (5.10) 

^ " n f(MHz) 

—8 ~7 

Thus, for 10 < n < 10 , the duration of the decay phase at 30 MHz is 

between 0.1s and Is, in excellent accord with the observed burst durations. 

As a further test, we note that because the decay duration t is less 
th 2 ui or of the order of the observed duration of elementary bursts, the 
source dimension is essentially limited to being less than or of the order 
of the distance traversed by a wave of group velocity = o^c in the observed 
duration. Taking “V 0.20, we find « 0.08 Rgif the observed duration 
is Is. This is an apparent size of 1.3 arc min, consistent with the 
apparent sizes of type III burst sources (Palmer and Lin 1973) . 

5. Polarization 

In the decay instability described in IV, the daughter 
TE wave has its electric vector along the electric vector of the pump 
wave, which in turn is approximately parallel to the magnetic field, which 
guides the beam. Thus the TE wave is emitted at large angles to the local 
magnetic field, as indicated by observations (see ^I) , and is linearly or 
highly elliptically polarized. 2\ny circular polarization observed 'in 
elementary bursts must thus be mostly impressed upon the waves during their 
propagation through the corona. Because all the striations in a single 
given elementary burst group originate in the same source volume, they 
follow the same iray paths and thus will exhibit the same sense and degree 
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of limiting polarization. Each elementary burst in a chain, however, 
originates in a different volume and so the limiting polarizations may vary 
among elementary burst groups, which traverse different ray paths. 

6. Energetics 

As discussed above, type Illb bursts are beeuned mainly at 
large angles to the local magnetic field. Therefore, knowing the 
source-region diameter D^, we can estimate the energy requirement by 
considering the power flow across the sides of a cylinder. The received 
power is 


P = W irDc. 6r (erg cm“^ s~^) , 

y g 


( 5 . 11 ) 


where r is the distance from the source to the observer, 6r is the radial 

dimension of the source region, and W is the energy density of TEM waves. 

For a single striation 6r = ^ L; using Newkirk's active streamer modfel 

f -3 

(fig. 9) , we have L - .5 Rg at 30 MHz, and fir =: 1.5 x 10 R^. Thus, 


0 


taking V = 0.1 c, r = 215 R_, and Dg = 0.08 R. in accord with our previous 
y 0 y 

estimates , we obtain 


-2 -1 

p 20 W erg cm s 


From the arguments above, where we calculated t, we see that W should be 
taken as the saturated energy density of the pump waves. Thus 

p =; 20 "e'^e cm“^ s ^ . 

Taking n = 10^ (f - 30 MHz) , = 2 x 10^ K, we get finally 

e e 


„ . T --2 4/3 -2 -1 

P 7 X 10 n erg cm s 


( 5 . 12 ) 


The received intensity of type III emission at decameter frequencies is 
typically of order 10~^^ erg cm“^ s“^ Hz“^, and in a bandwidth of 100 kHz 
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this becomes erg cm”^ s~^. Type Illb bursts are characteristically 

one to ten, times as bright as the following type III emission. Thus from 


(5.12), we require 


4/3 


> 


~9 

10 o. 


where 1 < a < 10. Taking 0.05 < 5 < 0.1 and a - 10, this becomes 

> 2 X io"ii -10"^°, 
or 

n > Ci-3) X lo**®, 


in excellent accord with the values of ri required by the other facets 
of the theory. 

7. Drift Rate 

In ^ I it was mentioned that in a very few cases , elementary 

bursts were observed to drift in frequency with rates |fj < 150 kHz s~^. 

the average rate being “ -18 kHz s ^ (de la NoS 1975) . Such frequency 

drifts are unusual and are not intrinsic to the simple theory we have formulated 

here. On the other hand, the fact that 99.8% of the obseirved drifts are 

negative indicates that when they occur, they are a systematic effect. 

One possibility is that such frequency drifts may be accou.'’ted for in a more 

detailed kinetic treatment of the theory that does not completely 

separate the time scales for the various instabilities. Another possibility 

which fits naturally into the context of our theory, is that the frequency 

drifts are due to variability in the source region itself. 

To see this, we note that for the daughter TE waves, the group velocity 

V = a c; from Figure 6 we can see that is independent of n — and therefore 
g T T ^ 

of either n or n, — and depends only on V = V /V, . From Eq. (5.2) we see 
e b e 9 

that the a^-dependence of a striation group depends on Therefore 
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variations in either the electron temperature or the beam velocity across 
the source region will lead to slight variations in the refractive 
indices — and therefore the frequencies and group velocities — of the 
daughter waves. (We note, however, that these variations are not sufficient 
to mix the different elements of a striation group.) 

Consider then the following simplified model of ray propagation in and 
near the source region. The refractive index goes to unity over a distance 
S which we expect to be roughly comparable to a density-gradient scale length L. 
We assume that all the waves propagate in one direction at near their 
initial group velocity, to the "cross-over” point S, after which they 
propagate at c. If S » D^, we may neglect the variation of position 

within the source region, and all the waves travel essentially the same 
distance. The time for a ray with frequency f and group velocity to reach 
S is given by t(f) = S/V . Then the frequency drift measured by an 
observer at (or beyond) S is 


dt 


V 


3 

T 


S dV /df 

g 


where we have used V = a_ c and 2irf = 03 (1 + a_^/2) . To estimate 

g T e ip / ' 


-1 


whether this is a reasonable model, we set f = 30 MHz, df/dt = -30 kHz s , 
so that 


a_3 5. ~ 10 

T S 


-3 -1 

s 


Then for 0.1 < a < 0.2, we find 3 x 10^^ < S < 2 x 10^^ cm, in reasonable 

f\j i\j ’ '\i '\j ' 

agreement with the assumption 'V/ 0.1 « S and S > L. The; few cases 

(~0i2%) in which df/dt is positive may be attributed to instances in which 


either is abnormally large or S abnormally small, so that 'v- S and 
variations of the path lengths to the cross-over region become important: 
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then source regions producing lower-group-velocity waves may be closer to 
the "cross-over” region. 

Finally, we may estimate the order of magnitude of the source 
fluctuations necessary to produce the observed drift rates according to 
this model. The indicated total variability in frequency across the 
source is given by 




where t is, as eibove, the burst duration. But 2irf = ‘‘ig (1 + a^^/2) , and so 

Af . AF 

— aAa ^5- 

X a§ 

Taking a - ? - 0.1, we have da^/dC - 2 (c.f. Figure 6), and AC C AT^/T^ if 
the variation is owing to temperature fluctuations, while AC CAU/U owing 
to variations in the beam velocity. In either case, fluctuations of 
< 1-10% across the source region may account for drift rates of the order 
observed. 

B. Chains of Elementary Bursts 

The above exposition has concentrated on the properties of 

elementary bursts. We now consider the reasons why many such burst groups 

may occur sequentially at discrete frequencies, more or less widely spaced, 

to form the chain that we designate as the type Illb burst. 

Again we return to inequality (1. 1) , which is the fundamental 

condition that must be satisfied in order that our mechanism may work. 

Above, in considering the frequency separation between striations in cin 

—8 —7 

elementary burst group, we found that density ratios of order 10 < r) <10 

were indicated by the observations. We now show that ratios of this order 
are in fact necessary in order that the linear instability discussed in 
be driven to saturation. 
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The linear instability may, of course, be initiated by beams 
with n « 10 i provided (1,1) is satisfied. Indeed, for collision 
frequencies 10“® 0 )^, we may have ^ for ry as small as 10”^^. 

Waves are unstable, however, at phase velocities in a narrow interval 
slightly below the instantaneous beam velocity U(t). In this velocity 
interval, the beam distribution function is either zero or rising to its 
peak at U. Therefore, in order that the unstable waves may be driven to 
saturation, we require that the evolution of the beam through this velocity 
interval take longer than the wave growth time injection 

at 3t=0, the time at X > 0 over which the beam evolves from a peak velocity 
U to a peak at D-AU is given by 


(At), - X 

beam — 
U 


(f)- 


. 1/3 

Setting Au = = p U, we then require 


(At) 


beam 


(At) 


sat 


V3 ,,, 

p x/li 

N-VryT/T'^ 

sat' e 


2/3 

^ ^ » 1 . 

N . U 
sat 


(5.13) 


8 “1 

Evaluating (5.13) at typical parameters of = 2 x 10 s (f - 30 mHz), 

-9 

X 0.6 R-, U = 0.5 c, and N = 10, we find p » 10 ■ . As soon as the 

beam distribution simultaneously satisfies inequalities (5,13) and (1.1 ) , 

however, the mechanism should begin to operate. Thus, allowing for variations 

in the coronal density model, inequality (5.13) is in accord with the values 
•“8 “7 

10 < p < 10 inferred earlier. 

'\t 'Xi 

Because the beam density builds up gradually at any point r, while 
the ambient coronal density decreases sharply with altitude, there is some 
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minimum radius at which inequalities (5.13) and (1.1) are first satisfied 
simultaneously. Of course, owing to variations in both the coronal density 
structure and the beam velocity profile, r^ may vary from event to event. 

Let us consider, however, what happens to the beam at r^ during the prpduction 
of the first elementary burst. 

The energy of the trapping wave, given by Eq. (2.14)., is derived 

at the expense of the beam particles. As a consequence, the particles 

1/3 

spread in velocity by an amount Av = q U (Drxmimond et al* 1970)^, In 
the single-wave theory, which does not consider the generation of sidebands, 
the beam and the wave exchange energy on the time scale so that the 

beam is sequentially reconstituted with a narrow velocity spread Vj^ 
and then spread again, the process repeating once every bounce time (O'Neil 
et al. 1971). With the growth of the sidebands, however, the beam particles 
are untrapped, because the beat frequency between the sidebands and the 
main wave is - u . This untrapping must occur over at least one bounce 

B 

period, and so after the particles are untrapped their thermal spread Vjj 

1/3 

is of order Av; i.e., Vj^ = U U. Then inequality Cl,l) is no longer 
satisfied; the beam distribution will be as schematically depicted in figupe 
11-b. 

Let us assume that the beam is- released from the trapping wave at 
the point r^ and at time t^; t = Q denotes the beginning of injection. 
Although it is not strictly necessary to do so, we assume that injection is 
still occurring at time t^» As long as injection is still occurring, the 
presence of particles of velocisty v at a point r means that particles with 
velocity v' > v are also at r, having been injected at times later than 
the injection time of the v-particles. On the other hand* at the head of 
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the beam the peak velocity U increases with increasing r. Consider the 
released beam (Figure 11) . Owing to the velocity dispersion in the beam, 
at a point r > r^ the particles at velocity U will again form the peak 
of the distribution function at the time t - (r-r^^j^j )/U, where r^^^j is the 
injection level. At r, however, the beam will again be narrower in velocity 


space than at r^, because the slower particles at r^ have not yet evolved 
to r and the faster particles present at (r,t) reflect the injected 
distribution, not the broadened one. We nov; ask: what is the minimum distance 
r , at which another elementary burst may be produced , and how many such 
bursts may occur in a chain? 

Assume for the moment that h (^j^) == n ; this assumption is 
examined below. Then the requirement (!•!) must be satisfied at r^ for 
a time of the order of t > the time required to saturate the sidebands. 


-5 

In Appendix A, we find 10 < t < 10 


-4 


Thus we require that the 


velocity-broadened particles at v < U not evolve to in the time t^j^this gives 


t < Ar n^^Vu, 
sb 


(5.14) 


where Ar = r,-r . Then the frequency separation (Af) between the 
1 o s 

elementary bursts at r^, r^ is approximately given by 

(Af), 


f (r ) 

o 



dm 


sb 

£i — 

e 

Ar > 

'V 


0) 

e 

dr 

2 1 rr 
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'u 10 , 

, U = 0.5c, 

o 

II 


(5.15) 


.-5 


.-4 


we find 


5 X lo”"* < 


(Af) 


'u f (r ) “V 
o 


® < 5 X lO”^ 


(5.16) 
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The range (5.16) is in general agreement vrith observations, 
although frequently the elementary bursts may occur so close together in 
frequency that the chain appears to be a continuous burst with frequency- 
modulated intensity. The parameters used to evaluate (5.16) , however, 

/ary widely and so (5.16) is to be taken merely as a general estimate. 

We turn now to the second question ; how many elementary bursts 
may occur in a chain? 

Consider first the radial density profile of the corona. In 
Figure 12 we illustrate the three density models mentioned above , for which 
the scale lengths L and intrinsic optical depths were calculated earlier. 
For the sake of definiteness, let us consider the intermediate model by 
Riddle (1974) which is given by 


ne(p) 



(5.17) 


where p = r/Rg is the normalized heliocentric distance. The decanieter 
range 20-80 mHz lies between p = 1.3 and p = 2.0. At p = 1.3 the terms 
in p~® and p”^® are in the ratio 3:1, respectively; at p =2.0 the term in 
dominates by a factor of 12. Thus we may represent the coronal density 
in this region by 



(5.18) 


where 6 a(p)< 10. 

The beam, however, is being guided by the magnetic field, and 
so at comparable stages in the beam evolution at r and r^^ we expect the beam 


density to be given by 
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nb(r) = xib(r^) 



( 5 . 19 ) 


Therefore, for a given beam velocity U, the density ratio of the beam to 
the ambient plasma v/ill be 


/X , \ ! \2 - a (r) 

n(r) = n(rQ) / (5.20) 

v^/ 

Therefore, we see that the natural evolution will tend to make ri/ 
corresponding to a given value of U, ^ arger the greater the distance from 
the injection point. On the other hand, the spreading of the beam by the 
trapping wave at r^ must remove some fraction of the particles from the 
velocity range v = U; let us denote the relative magnitude of the 
distribution functions before and after spreading by 

|P,(U)| after spreading 

XE — . 

I (U) [before spreading 

Obviously X < 1, but we expect it to be of order unity. Then assuming 
an elementary burst to have been produced at Pq, Eq. (5.20) wilJr be modified 
at p > Pq to 

n(p) = Xn(Po) /p_'f • <5*21) 

Vo) 

As X'fe have seen above, in order for an elementary burst to 
occur, the inequalities p » lO"^ and V^/U < must simultaneously 

be satisfied. An estimate of the bound on the number of elementary bursts, 

B, given by the latter of these conditions may be inferred from (5.14) 
or (5.16), from which we may estimate B to be of order 1-10 over the decameter 
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range. Figure 13, from de la Noe and Boischot (1972), shows the observed 
distribution of B. It can be seen that most bursts indeed have B < 10. 

An estimate based on (5.21) may be obtained by assuming p =2.0 
ajid a ~ 1, over this distance range. Then, since the beam density will 
be modified by approximately the factor A in each burst, and since we 
require n (p=2) = r|(p=1.3), we may estimate B from 

^ x®(1.5)^ = 1, (5.22) 

Po 

from which B = -2/lnA. In Table 1 we show some values of B vs . A according 
to (5.21). Although it is impossible to estimate A precisely without making 
detailed calculations of the beam dynamics and evolution, we see that 
the range of observed values of B may be accounted for by quite plausible 
values of A . 

We stress here that we have made these estimates for a particular 
density model. The range of parameters associated with the various possible 
density models, together with the general variability of the parameters 
governing both the beam evolution and the instabilities, allow for 
correspondingly great variability among individual type Illb bursts. 
Nevertheless , we see that it is indeed possible to combine the general 
features of the beam evolution and the theory for the elementary bursts, 
to explain how the burst chain may be produced. 

Finally, we note that in the context of Eqs. (5. 14) and (5.21), the 
simultaneous requirements on n and Vj^/U may become difficult to satisfy 
simultaneously after many bursts have occurred or when the head of the 
beam has propagated to a region where the coronal density variation is 
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sufficiently slow so that ri (r) decreases, owing to velocity spreading. 
These considerations indicate why type Illb bursts are not observed below 
about 10 MHz. 
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Appendix A - Injection and Space-Time Evolution of the Beam 

Most type Illb bursts occur in conjunction with decametric continuum 
or dui:ing storms of type III bursts at decameter wavelengths, which are 
associated with noise storms at meter wavelengths. This association indicates 
that types III and Illb bursts in such storms occur when previously- 
accelerated particles escape from magnetically-confined regions to open 
field lines, and furthermore that this injection must be rather sharply 
delimited in time. The exact injection mechanism, however, is unknown. 

Let us then consider the following phenomenological model for the injection 
and evolution of the beam in the lower corona. The model, of course, is 
neither unique nor compelling, but merely provides one of several possible 
contexts in which the assertion of inequality (1-1) ivS plausible. 

We assume injection to take place at x=0, where the beam distribution 
function is taken to be of the form 

Fj^(x=0,v,t) = fo(v) I(v,t). (Al) 

We do not attempt to model an injection process that forms a Lorentzian 
distribution such as was assumed in^II, because that shape was used there 
only for its analytic convenience; as noted in^II, similar results 
concerning the beam-plasma instability may be found within a wide class of 
distributions, provided inequality (1.1) is satisfied. 

It is reasonable to consider fo(v) to be a monotonically decreasing 
function of v. The function I(v,t), which we shall call the "injection 
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profile," must rise over a time much shorter than the typical duration 
between individual bursts in a storm, which in general may be of the order 
of seconds. We take for I the form (Figure 14) 



0, t > T (V) . 

As we shall see, the value of t (v) is not terribly important for 
consideration of the early evolution of at x > 0, nor is the precise 
form of I for times t > t ( v) , provided I(v,t) reaches a maximum at t = x (v) . 

A particularly simple class of velocity -dependent injection profiles 
may be obtained by choosing 


x(v) = Tq 3 
a (v) = Oq g“^. 


(A3) 


where 3 h v/c and x^, a^, a are constants. In particular, we note that 
I becomes velocity- independent for the choice a = 0. In the limit 

->■ 0, I approaches the unit step fxinction H(t), and the injection is 
impulsive for all velocities. 

Using the forxns (A3), it is clear that if Ff, (v) is monotonically 
decreasing, F^^(x,v,t) will be sharply peaked near the value U(x,t) = 
3 } 3 (x,t)c such that 
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(A4) 


Owing to the modulated gaussian envelope of the rise of I(3,t) and to 
the increase of a (3) as 3 decreases, the thermal width Vj^ of the beam 
satisfies 


^(x,t) <1 3j^(x,t) - 3j^(x,t - Qq3^) 
c 

= Oq 3b^'^^ (x,t) 

r 

t 

where we have taken t = x/3c = tj,/3. Then 


V. 


_b 

U 


3b 


a+1 


and inequality l-l, becomes 



(A5) 


(A6) 


For n - 10~® - 10“"^, this condition is fairly easily satisfied. We note 
that in the decameter range, t - Is, and must be considerably less 
than 1 s because the overall time scale of the injection is also < 1 s. 
Thus for the sake of illustration, let us take (Uq/To) - 0.1 and a = 2. 

Then for n = 10'’^, wc would require 3 ^ 0. 3, quite epnaintent with our 
earlier considerations in Moreover, we may see that a wide range of 

plausible values for the parameters may be chosen to satisfy (A6) . 
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Finally, we note that because we have used the one-dimensional form- 
ulation of the beam-plasma interaction (c.f. ^ II), inequality (1.1) actually 
applies to the projection of the distribution function on the direction of 
the magnetic field. Therefore, even if the beam were strictly mono- 
energetic, the spread in pitch angles A0 would be required to satisfy 


AG < n 


(A7) 


in analogy with (1. 1) .Therefore, we must inquire whether pitch-angle 
scattering will broaden the beam in a time shorter than that required for 
the sidebands to grow. The mean-square pitch-angle scattering of an 
electron due to Coulomb collisions, after traversing a distance L in 
the plasma, is given by 


<(A0)^> 


8-?r He e'^ 
m^ v^ 


InA 


L, 


(A8) 


O *3 

where v is the electron velocity and A = (3/2 e"^) (T^ /it n^) . 

Evaluating (A8) at = 30 MHz (n^ - 10^) and v = 10 - .2c, we find 

12 

that a path length L = 5 x 10 cm is necessary for a mean root-square 
scattering of 10“^ rad. Thus, providing the projected distribution function 
satisfies inequality (1.1) at any time. Coulomb collisions are absolutely 
negligible as a mechanism for invalidating the theory. 

The other possible source of pitch-angle scattering is magnetic 
fluctuations. The effects of such scattering are somewhat difficult to 
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determine accurately. In the quasi-linear theory (Jokipii 1971) , the 
pitch-angle scattering coefficient is given by 


D, 


W 


AMh. = (Ae)2 


(A9) 


where e = Lg/Rg, ji = cos 0, and A = (5 b/B)j^s is the rms value of the 
field fluctuations. Lg .is the correlation length of the field fluctuations 
Rg the electron gyroradius, and At is given in correlation times: 


At = vt/Lg. 

A precise assessment of the scattering coefficient is difficult because 
the parameters in (A9) are, of course, not well known very near the sun. 

Let us assume that the scale length Lg is '^^10^ cm, the scale length of 
supergranulation, and that the field strength is of order 1 gauss. Then 
for V .3c, e'vlO^, and (A9) becomes independent of e. Writing 
<(Ay)^> - (1-y^) < {A6) ^>, taking |y|;vlf setting A6 we find 

the scattering time Tg^^tt 


^scatt 



(AlO) 


Setting A ''^0.1 (the approximate upper limit for validity of the quasi- 

-8 -4 

linear theory), v =: .3c, and y == 3 x 10 , we find - 10 s. 

Assuming ten e-foldings are necessary for saturation of both the trapping 
waves and the sidebands, we find tgj-Qy^th ~ 2O/y^/^03g ~ 3 x lO”^ - lO”'^ s. 

Thus, within the context of this rather speculative calculation, pitch-angle 
scattering by magnetic irregularities may be marginally sufficient to compete 
with the tendency of the beam evolution to satisfy (1.1). 
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Appendix B r. REMARKS ABOUT THE SINGLE-WAVE THEORY Al'ID SIDEBAND INSTABILITY 

The theoiry we have presented above invokes three distinct processes 
that occur sequentially on well-separated time scales: the strong beam- 

plasma interaction, which we have described in the context of the single- 
wave theory; the parametric sideband instability; and the parametric 
decay instability. Although each of these elements of the theory has 
been discussed more or less extensively in the literature, there remains 
some degree of contradiction between theoretical and experimental results, 
and our understanding of the individual processes is not yet complete. 
These remarks apply particularly to the sideband instability. Further- 
more, to our knowledge, the sequential evolution of the three processes 
such as we have described above has never been observed in laboratory 
plasmas. Therefore, although at the least we may consider the theory 
merely as a model which successfully organizes the observed phenomenology 
of the type Tllb bursts, it seems desirable to discuss briefly the 
present state of evidence bearing upon the elements of the theory. 

The single-wave theory was tested in a comprehensive experiment 
by Gentle and Lohr C1973) , who verified the theory in all respects up 
to the time the main wave passes through one amplitude oscillation. 

At that time, they observed the main wave to dissipate into an incoherent 
turbulent spectriam, while the particle distribution function became very 
broad. It is implicit in our mechanism, of course, that the main wave 
remain coherent long enough to amplify sidebands to a significant level. 
Assuming that the sidebands grow to levels comparable to that of the 
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main wave (not necessarily a good assumption, as we shall discuss 
below) , estimates in Appendix C indicate that of the order of ten 
e-folding times Ygj^ ^ required. 

The reasons for the breakup of the trapping wave in the experiment 
of Gentle auid Lohr are not well understood (K.W. Gentle, private 
communication) , and different effects may dominate in different parameter 
regimes. Some possible mechanisms that may play a role, for example, 
are the oscillating two-stream instability and three-wave coupling. 

The experiment may also be influenced by gradients. In addition, it 
is unclear how the experimental results may scale with ri. The breakup 
of the wave is probably not due to demodulation by the more slowly- 
growing waves in the initially unstable wave packet, however, for it 
is found both theoretically and experimentally that the saturation of 
the fastest-growing wave suppresses the further g' owth of nearby waves 
(DeHeef, Malmberg, and O'Neil 1973). 

Gentle and Lohr did not attempt to observe sidebands. Beam-plasma 
experiments in which sidebands were observed were conducted by van Wakeren 
and Hopman (1972), Mizuno and Tanaka (1972) Jungwirth, Piffl, and 
Ullschmied (1974), and Nyack and Christiansen (1974). Thus, there is 
direct experimental evidence that sidebands i amplified in beam- 

plasma interactions, the initial stages of which are described by the 
single-wave theory. 

With regard to the sideband instability itself, however, there is 
general disagreement between experiments and theory and among different 
experiments. Existing theories may be classed as either parametric 
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(Kruer et al. 1969; Goldman and Berk 1971; Wong 1972) or "quasilinear" 
(Bud'ko, Karpman, and Shklyar 3.971; Brinca 1972). In the parametric 
case, only the interaction between resonant electrons and a single 
driven mode is considered, while quasilinear theories consider the 
interaction with other modes during the evolution of the distribution 
function. There are also two general types of experiments — those 
in which a beam amplifies the trapping wave, and those in which a 
large-amplitude wave is launchea into the plasma. 

Parametric theories of the instability predict growth of both 
upper and lower sidebands. In the quasilinear theory of Bud'ko et al. , 
the lower sideband (with slightly faster phase velocity) dominates, 
while Brinca finds both upper and lower sidebands, although he argues 
that the faster sideband should dominate, owing to Landau damping of the 
slower one. Ejq>eriments , on the other hand, generally observe either 
one or the other sideband, van Wakeren and Hopman, Mizuno and 
Tanaka, and Jungwirth et al . all observe the lower sideband, while 
Nyack and Christiansen observe the upper one. These latter authors, 
however, argue that these results are consistent if the slower sideband 
is the unstab] e one, because in the experiments of van Wakeren and Hopman 
and Mizuno and Tanaka the linearly unstable wave is near the upper 
hybrid frequency and the upper sideband corresponds to the slow space- 
charge wave, CThis comment appears to apply also to the experiment of 
Jungwirth et al.) The above experiments all involve beam-plasma inter- 
action. In the launched-wave experiments of De Neef (1974) and Van Hoven 
and Jalms (1975) , the lower (faster) sideband was observed. In both of 
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these experiments, however, coupling between sidebands is also observed. 
Jahns and Van HOven (1975) attribute this coupling to a passive four- 
wave mode coupling process, thereby arguing against the parametric 
interpretation of the instability. Thus, whether or not the sideband 
instability is parametric, experiments indicate the possibility of 
amplification of both upper and lower sidebands. 

A further discrepancy between theory and experiments arises when 
we consider the amplitudes of the main wave and the sidebands during 
the instability. In we assumed that the sidebands grow to amplitudes 
comparable to that of the main-wave saturation amplitude. This 
assumption does not necessarily violate conservation of energy, because 
the energy in the sidebands may come from the particles. In addition, 
all the theoretical treatments we have cited assume the main wave remains 
at constant amplitude. (Kruer and Dawson D-970;3 , though, did computer 
simulations in which the energy of the main wave was depleted as the 
sidebands grew.) In both theory and experiment, however, the amplitude 
of the main wave oscillates over a factor of about five during the bounce 
period, and corresponding oscillations are observed for the sidebands. 
Nevertheless, this fact is not crucial to our theory, because we considered 
the wave amplitudes in the context of the energetic requirements of type 
mb bursts, and found in that for the most demanding case (a - 10) 
these requirements were well met. In the same vein, we note that the 
experimental results of Van Hoven and Jahns are that the amplitude of the 
main wave decreases about two orders of magnitude while the sideband 
grows, both waves asymptotically approaching the same amplitude. Jungwirth 
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et al., also found the sidebands to approach the amplitude of the 
main wave. Again, there appears to be no essential difficulty for 
the theory here, since our estimate of the energy requirements was 
shown to be well satisfied. 

Finally, we note that the growth rate for the sideband instability 
that we calculated in III is Ygj^ ~ a value which we feel pushes 

against the limits of validity of the linear theory of Goldman and Berk. 
The experimental results of Van Ho-o-ens and Jahn and of Mizuno and 
Takanaka, however, indicate > Wg. 

Taken in the aggregate, these remarks show the wide divergences 
between theoretical and experimental results concerning the sideband 
instability. Part of the discrepancy doubtlessly lies in the fact that 
the theoretical treatments are for the case of temporal instability, 
while the observations are of spatial instability and include the effects 
of inhomogeneities, magnetic fields, and plasma boundaries. Nevertheless, 
it seems fair to assert that although the validity of the single-wave 
theory and the existence of the sideband instability are beyond doubt, 
a comprehensive theory of the latter is yet to be attained. In our 
theory of type Illb bursts we have used the sideband instability theory 
of Goldman and Berk, in part because it is one of the simplest and most 
tractable theories extant. In addition, of course, the plasma 
conditions of the solar corona are much closer to the ideal conditions 
assumed in theoretical treatments than to the conditions of laboratory 
experiments, in that the corona is essentially uniform and infinite on 
the scales of the relevant interactions. 
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In this regard, we note that for the time scale of the type Illb 
burst, which in we showed is detemnined by the decay rate y^, an 
experimental system such as that of Gentle and Lohr is of order ten 
thousand times too small to permit the sequential occurrence of the beam- 
plasma interaction, sideband instability, and parametric decay — the 
three fundamental elements of the type Illb burst. For this reason, 
we do not expect to see phenomena such as the type Illb burst in the 
Iciboratory; the solar corona is a much more ideal plasma "laboratory" 
than can be created artificially. 

Finally, we note that the qualitative details of our theory are 
relatively insensitive to the exact theory employed to describe the 
sideband instability, for three reasons. First, all theories predict 
the sidebands to be separated from the main wave by approximately the 
bounce frequency, and this was the essential quantitative result used in 
^V. Second, the time scale for the type Illb phenomena is essentially 
determined by the time scale for the parametric decay instability, 
because this is the longest time scale in the inLeraction regardless 
of the precise growth rate of the sideband instability. Third and 
finally, the energetic requirements of the type Illb bursts appear to be 
easily met by the theory notwithstanding the assumptions we have made 
regarding the amplitudes of the main and sideband waves. 
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Appendix C - THRESHOLD CONDITIONS FOR THE PARAMETRIC DECAY AND UP-CONVERSION 
INSTABILITIES 

In ^IV we stated, without proof, that the only parametric instability 
for which the threshold pump intensity is exceeded in our theor^y is the 
parametric decay instability described by Lashmore-Davies (1974a, b) . 

In this APPsndix we support this assertion by estimating the enhanced 
TE-wave damping rate in the source region, and computing the ratios 
of to for both 5 < 0 and 6 > 0 ^c.f. Eg. (4.-4)J . 

First, we note that the electric field of the TE. wave is along the 

I -A. 

direction of the pump wave field, so that k.^, As we saw in ^V, the 

magnitude k^ = arpWg/c 5 0 . 2 , But kj^ == Ug/U, so that k,j,/k^ < 0.2 U/c « 1. 
Therefore, according to the resonance conditions (4.-1) » ~ ~ ^L‘ 

allows us to write 

== ^s ~ ^ , (Cl) 

“e u ^ 

where y = mg/mp is the mass ratio of electrons to protons, and Cg = y’^V^ 
is the ion sound speed. In an isothermal plasma (T^, = T_) , the ion sound 

e p 

wave is heavily damped; ignoring factors of order unity, we may take 
= Wg for the purpose of making estimates. 

The damping rate Yiji the background plasma is initially given by 
Y^ ~ J^e^e ^s the niimber of particles in a Debye sphere 

[Xg 5 (T^/4 tt n^e^) However, as the decay progresses, the buildup 

of ion sound turbulence leads to enhanced resistivity of the plasma to 
waves near the plasma frequency (Dawson and Oberman 1963; Dawson 1968). 

Thus, the threshold changes in time, and we should use the maximum 
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value of Y attainable. For an energy density W of sound waves, in 
i s 

which the typical wavenumber is k , the damping rate ytp n^^y be estimated 
(Papadopoulos 1974) as 




(Wg/ne^) “e 
k 2 X^2 


(C2) 


Tlierefore, we need to estimate the maximum level of W^. 

Denote H W^^/n^Tg, S E Wg/ngT^, where is the energy density 
of the main wave. In the wavelength interval k^ - < k < k^^ + 6k, the 

initial noise level Rq( 0) is given by 


RqCO) = _2 (koXe)3 ^ ^ 2 , (C3) 

% ko Np 

where to get the RHS we have used Eq. (2-10) . Therefore, using Eq. (C3) 
and Y^(0) - t 0 e/ 2 Nj^, we have 

K(0) = ^1/3 ^ ;j^Q ^2 ^1/3 

K^,(0) 4 

Using Eq. (2.14), the saturation level is 



sat 




(C5) 


Denote by x the total time the main wave takes to grow from the initial 
noise level to saturation. Further, we denote t)y duration between 

the time the wave reaches K (0) and saturation. Then x. is given by 

C X 

Ksat = = 5 = exp (2^ x^) , (C6) 

Kc(0) 
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while T is defined, using Eqs.(c4) and (C6) , by 


Ksat = = expC2-^x) 


(C7) 


K, 


2C' 


Denoting the decay rate found from Eq. (4.2) by the sound level 

is then given by 

®sat "" ^ exp{2y^T3_) , 

where S (0) - R{0) . Then 


\,sat = y,p(0) exp(2y^TQ^) = 2K(0) ^xp (2 


(C8) 




where is the beam-plasma growth rate found in^II, For 

an overestimate, we take y,^ - then Eq. (C8) becomes 



= 2K 


w, 


5^ 

L 


4/3 


y*»€/nV3 


(C9) 


where we have used Eqs. (C4) and (C6) ■ Substituting Eq. (C9) into the 
first of Eqs (4.4) to determine the decay threshold at saturation, we 
find 


^sat 
^c , sat 


Npn 

liy^ 


5Npn 


4/31 


1/3 

-y C/n ^ 


(CIO) 


^sat' given by (CIO) are plotted in Figure 4. It is 

seen that under conditions of interest to our theory, K greatly 

s d.'C 

exceeds the decay threshold but is in general much smaller than the 

up-conversion threshold K = 

sat 


2 . 
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TABLE 1 

Niomber of possible elementary bursts B for various values of 
X: B = -2/S,n X. 


X 

B 

0.50 

3 

0.60 

4 

0.70 

5 

0.80 

9 

0.85 

12 

0.90 

19 

0.95 

39 
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Figure 1. 


Figure 2. 


Figure 3, 


Figure 4. 


Figure 5. 


FIGURE CAPTIONS 

Dynamic spectra of type Illb bursts. (a) A type Illb burst 
followed by a faint type III burst on August 8, 1970. (b) Another 

type mb burst followed by a type III burst on August 12, 1970. 
Vertical dotted lines are lightnings. Horizontal white lines are 
interferences . 

Two-dimensional phase space in the frame moving with the 
phase velocity of the trapping wave. 

(a) Maximvun growth rate and corresponding wavenumber 
shift and real frequency shift for the sideband dispersion 
relation of Goldman and Berk (1971) , plotted against ^ 
for £ = 2. All quantities are normalized to 03^. 

(b) Asymptotic C^-^-O) values of the maximum sideband growth 
rate and corresponding frequency shift, normalized to 

as a fiuiction of ri. 

o> 

(a) Saturation amplitude K_ of the main wave as a function 
of V, for various ri. 

'\i f\j 

(b) Ratio of K .to the parametric decay threshold K 

S C f S3.W 

Maximum growth rate of the parametric decay instability, as 
a function of the phase velocity of the trapping wave, 

= ?. The equivalent decay rates of the sidebands may be 
obtained by multiplying by (V+/V^)^. Dashed curves--for a 
monochromatic pump. Solid curves — corrected for finite frequency 


spread in the pump. 
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Figure 6. Index of refraction of the TE' waves resulting from the 
decay of longitudinal waves at the effective normalized 
phase velocity V. This curve is independent of n. 

Figure 7. Bandwidth of the decay wave spectrum as a function of V, 
for various n. 

Figure 8. Effective phase velocities V+ of the sidebands, as a function 
of Vq, for various ri* 

Figure 9. Scale length L of the coronal density gradient for various 
density models, and the corresponding intrinsic optical 
depths T^. 

Figure 10. Optical depths of the decay daughter TE' waves, as a function 
of V^, at 30 MHz and 100 MHz. Tq,t~, and refer to the 
main wave, red (lower), and blue (upper) sidebands, respectively. 

Figure 11. Schematic illustration of the beam distribution function 
Fj^(v): (a) before trapping; (b) after being untrapped. 

Figure 12. Coronal density models of Leblanc et al. (1973), Riddle (1974), 
and Newkirk (1967) . 

Figure 13. occurrence histogram of the number N(B) of type Illb bursts 
vs the number B of elementary bursts. 


Figure 14. 


Schematio illustration of the injection time-modulation 
profile l(v,t) . 
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